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Abstract 

Optical vortices arise as phase singularities of the light fields and are of central 
interest in modern optical physics. In this paper, some existence theorems are estab- 
lished for stationary vortex wave solutions of a general class of nonlinear Schrodinger 
equations. There are two types of results. The first type concerns the existence of 
positive-radial-profile solutions which are obtained through a constrained minimiza- 
tion approach. The second type addresses the existence of saddle-point solutions 
through a mountain-pass-theorem or min-max method so that the wave propaga- 
tion constant may be arbitrarily prescribed in an open interval. Furthermore some 
explicit estimates for the lower bound and sign of the wave propagation constant 
with respect to the light beam power and vortex winding number are also derived 
for the first type solutions. 



1 Introduction 

Vortices have important applications in many areas of modern physics including con- 
densed matter systems, particle interactions, cosmology, and superfmids. Research on 
vortices in optics also has a long history and was initiated in as early as 1964 by Chiao, 
Garmire, and Townes [5] who explored some conditions under which a light beam can 
produce its own waveguide and propagate without spreading. They described such phe- 
nomenon as self-trapping, attributed it to light propagation in materials whose dielectric 
coefficient increases with field intensity in the context of high-intensity light beams such 
as lasers, predicted marked optical and physical effects, and suggested the occurrence 
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of optical vortices. Such vortices have since then been observed in numerous studies 
[2, 3, 27, 29, 32] and become a much pursued subject in optics [6, 9, 11, 19, 22, 31, 33] 
both theoretically and experimentally. (For a vast literature up to 2005 and for a descrip- 
tion of the rich features and profound applications of optical vortices, see the nice survey 
article by Desyatnikov, Kivshar, and Torner [10]. See also [16] for a more recent survey 
of the subject in a broader perspective.) As waves, light propagation may be described 
by a wave function. At certain spots of space, the intensity of the waves vanishes and 
the phase of the waves cannot be defined. Thus, such spots are phase singularities which 
were recognized in the comprehensive work of Nye and Berry [23] as crucial character- 
istics of general wave motions in which vortices are present. These phase singularities, 
also referred to as dislocations or defects of waves, are the centers of vortices, around 
which energy concentrates. In the context of light waves, vortices are centered around 
vortexlines and light waves are twisted around the vortexlines. The twisting arises from 
the phase ambiguity around a vortexline and is of a topological nature. The twisting 
centers are exactly the vortex cores where light waves cancel out leading to darkness so 
that light intensity measured at any cross section vertical to a light beam axis should 
display concentric ring-like patterns around the dark core. Light beams of such structure 
are also vividly termed "helical beams" [4]. In optics research, a fundamental prototype 
situation is when the light waves are described by a complex-valued wave function gov- 
erned by nonlinear Schrodinger equations [1, 7, 17, 18, 20, 21, 26, 28]. These theoretical 
studies provide a broad range of interesting analytic problems related to the existence 
and properties of optical vortices for mathematical investigation. 

Our aim in the present work is to obtain some existence theorems for the optical vortex 
solutions explored by Salgueiro and Kivshar in [28] through a study of the normalized 
nonlinear Schrodinger equation 

i^ + V 2 ^+(V + s\^ = 0, (1.1) 

where ip is a complex-valued optical field propagating in the (longitudinal) z-direction, 
is the Laplace operator over the (transverse) plane of coordinates (x,y) which is 
perpendicular to the z-axis, V is an external potential function (cf. [18, 30]), and s = ±1 
is the sign symbol indicating either a focusing or defocusing situation [17] which is taken 
to be +1 (focusing) in [28] and will be our main focus (the defocusing case s = — 1 will be 
seen to be straightforward). The interest of (1.1) actually goes beyond nonlinear optics. 
For example, it also arises in the study of the Bose-Einstein condensates [1, 11, 15, 30] 
and is referred to as the Gross-Pitaevskii equation. An important simplified situation 
that allows optical vortices to present is when V depends on the radial variable only, 
V = V(r),r = \/ x 2 + y 2 . In this situation one may expect to find an n-vortex solution 
of (1.1) of the form [28] 

i> = ip(r, 9, z) = u(r)e iind+(3z) , (1.2) 

where r,8 are polar coordinates over M 2 , u(r) is the radial profile function which gives 
rise to the intensity of light waves, integer n e Z is the winding number or vortex charge 
of the vortex solution, and f3 e M is the wave propagation constant [28]. This ansatz 
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describes a vortex wave centered around the z-axis where r = and propagating along 
the z-axis. Inserting (1.2) into (1.1), we arrive at the following n- vortex equation 



of cubic nonlinearity. The presence of the vortex core at r = requires u(0) = 0. As in [28] 
(for n — 1), we are interested in ring-shaped vortices so that light intensity concentrates 
around the vortex core which suggests that u(r) may be assumed to vanish at a sufficiently 
large distance. Mathematically, this indicates that we may impose another 'boundary' 
condition, say u(R) = 0, at a certain distance R > away from the core of the vortex 
as seen in the numerical results of the work [28]. Thus, the problem of the existence of 
optical vortices is reduced into a two-point boundary value problem with undetermined 
parameter /3 and prescribed R, for any given n G Z. To tackle this problem, we shall use 
the methods of calculus of variations. Our methods allow us to obtain two types of results 
for the nontrivial focusing case s = +1. The first type of results rely on a constrained 
minimization approach. The nature of minimization leads us to obtaining positive-valued 
solutions in the open interval (0, R) and that the propagation constant j3 arises as a 
Lagrange multiplier due to the constraint. The second type of results are obtained from 
searching for saddle points of the action functional associated to the problem. We will see 
that, in this latter case, there is no assurance that the solutions must stay positive-valued 
but the propagation constant (3 arises as a prescribed quantity. 

In the next two sections, we shall concentrate on the focusing case when s — +1. 
In Section 2, we formulate the problem of existence of optical vortices as a constrained 
minimization problem, state the main existence results regarding positive solutions, and 
then present the proofs. We will see that the propagation constant j3 arises as a Lagrange 
multiplier which is ensured to be negative when the vortex charge n is sufficiently large. 
We will also derive some lower estimate for (3. In Section 3, we treat j3 as a prescribed 
quantity and use a mountain-pass theorem approach to establish the existence of solutions 
for any R and vortex charge n. In particular, we show that the propagation constant j3 
may assume any prescribed value in an explicitly given interval. In Section 4, we briefly 
discuss the defocusing case when s = — 1. 

2 Vortices via constrained minimization 

As described in the previous section, we will be interested in 'ring vortices' such that 
(1.3) is considered over a bounded interval (0, R) (R > 0) so that u vanishes at the 
two endpoints of the interval. As mentioned earlier, we will mostly concentrate on the 
nontrivial case, s = +1. Thus, our problem is a two-point boundary value problem 




(1.3) 




(2.1) 
(2.2) 



u(0) = u(R) = 0, 



for which the parameter (3 G R arises as an eigenvalue of the problem. 
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In order to approach the problem consisting of (2.1) and (2.2), we write down the 
action functional 



I(u) = \ f R iru 2 r + — u 2 -rV{r)u 2 - T -u 4 [> dr. (2.3) 

2 .In \ T 2 " 



and the constraint functional 

P(u) = J \ip\ 2 rdrd9 = 2tt^ ru 2 dr, (2.4) 

which measures the beam power [28] of the vortex wave. Thus, to get a solution of 
(2.1)-(2.2), it suffices to prove the existence of a solution to the following constrained 
minimization problem 

min {I{u) \ueA, P{u) = P } , P > 0, (2.5) 
where the admissible class A is defined by 

A = {u(r) is absolutely continuous over [0, R], u(0) = u(R) = 0, E(u) < oo} , (2.6) 
with 

E(u) = \J Q |™r + \u 2 + T -u^ dr, (2.7) 

being the 'energy' functional, Pq is a prescribed value for the beam power, and /3 arises 
as the Lagrange multiplier. Note that the finite-energy condition E(u) < oo is imposed 
only to ensure that all the terms in the indefinite action functional (2.3) stay finite. 

For convenience, for a function / of the variable r, we interchangeably use f r and /' 
to denote its derivative. We will also need the following decomposition and notation 



V = V+-V-, V ± = max{±V, 0}, 

= max{y ± (r) |r e [0,R]}, 

V = max{|F(r)| | r e [0, R}} = max{V + , V^}. j 
The main results of this section may be stated as follows. 



(2.8) 



Theorem 2.1. For any nonzero integer n and a given continuous potential function 
V(r) over the interval [0,R] (R > 0), consider the two-point boundary value problem 
(2.1)-(2.2) governing an n-vortex wave solution of the nonlinear Schrddinger equation 
(1-1), propagating along the z-axis with a propagation constant (3. 

(i) The problem always has a solution pair (u, 0) with u(r) > 0, r G (0, R), and /? 6 i, 
so that the associated beam power enjoys the bound P{u) < Air\n\. In fact, such a solution 
may be obtained through solving the constrained minimization problem (2.5) assuming 
Pn < An\n\, from which (3 arises as a Lagrange multiplier. 
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(ii) Let (u, pi) be the solution pair of the problem (2.1)-(2.2) obtained in part (i). Then 
P has the lower bound 

7 P 12 

^-^~^-^ (1 + n2[21n2 - 1]) - (2 ' 9) 

(Hi) Let (u,P) be the solution pair of the problem (2.1)-(2.2) obtained in part (i). 
Then P < if \n\ is sufficiently large so that 



\n 



(pi -)h 

> I ^ + max{rV+(r) \r G [0,R]} \ . (2.10) 



(iv) The problem (2.1)-(2.2) has no nontrivial small-beam-power solution satisfying 
P( u ) < \ if the condition 



n 



>r\V + (r)-P), re[0,R], (2.11) 



is fulfilled. So, roughly speaking, the problem has no nontrivial small-power P and small- 
propagation- constant (i.e., \P\ is sufficiently small) solution over a small interval [0,R\. 

We now establish these results. 

(i) For any function u satisfying u(0) = 0, the Schwartz inequality implies that 

u\r) = jT 2u(p)u p (p) dp < 2 QT pu%p) dp) 2 QT ^ dp) * . (2.12) 
Thus, multiplying (2.12) by ru 2 , integrating, and using P{u) = P , we have 



Inserting (2.13) into (2.3), we obtain 

'<»> * 5 - 5) jf "^ + W-ik) £ T Ar -l p ° v °- < 2 -") 

In order to be able to find a suitable e > such that in (2.14) we have 



simultaneously, it suffices to assume that P satisfies the condition 

P <47r|n|. (2.16) 
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In this situation, we can find two positive constants Ci,C 2 , depending on e,n,P but 
independent of u, such that 

r R r R v 2 l 

/(u)>Ci/ ru 2 r dr + C 2 — dr-—P V . (2.17) 

JO JO r 47T 

Assume (2.16) and let {u m } be a minimizing sequence of (2.5). Then the coercive 
inequality (2.17) gives us the bound 

ri[u m ] r ) 2 dr+ / -u 2 m dr<C, (2.18) 
Jo r 

where C > is a constant independent of to. 

Since both functionals I and P are even, we have I(u m ) > I(\u m \) and P(u m ) = 
P{\um\), where we have also used the basic fact [13] that for any function u its distribu- 
tional derivative must satisfy \\u\ r \ < \u r \. In other words, we see that the sequence {u m } 
may be modified so that each u m is nonnegative, u m > 0. Thus we may assume that the se- 
quence {u m } consists of nonnegative- valued functions. It is clear that these functions may 
be viewed as radially symmetric functions over the disk Dr = {(x, y) G M 2 | x 2 + y 2 < R 2 } 
which all vanish on 8Dr. Moreover, since (2.18) holds, we see immediately that {u m } is 
bounded under the radially symmetrically reduced norm | | where 

•R rR 



\u\ 



/ ru 2 dr+ / ru 2 r dr, (2.19) 
Jo Jo 



for the standard Sobolev space W ' (Dr) since 

-r r R 



r u f 1 

/ ru 2 dr <R 2 -u 2 dr. (2.20) 
Jo Jo r 



Hence we may assume without loss of generality that {u m } converges weakly to an element 
u G WI' 2 (Dr) as m — > oo. Applying the compact embedding W 1,2 (Dji) — > L p (Dr) 
(p > 1), we see that u m — > u strongly in L p (Dr) as m — > oo. Of course, u is radially 
symmetric as well. Thus we may write it as u — u(r) which satisfies u(R) = 0. Moreover, 
it is clear that for any e G (0,R), {u m } is a bounded sequence in the space W l ' 2 (e ) R). 
Thus, using the compact embedding W 1,2 (e, R) — > C[e, R], we see that u m — > u as m — > oo 
uniformly over [e,R\. Besides, similar to (2.12), we have for any pair ri,r 2 G (0,i?), 
f'\ < ?"2, the inequality 



\ h / r r '2 ii 2 \ ^ 



\<{r 2 )-u'Jri)\ < 21/ r([u m ] r ydr) [I ^dr 



n / \J ri 

1 



< 2<?5 



i(J"^rJ', (2.2!) 



where the constant C > is as given in (2.18). Letting m — > oo in (2.21), we arrive at 



\u 2 (r 2 ) - u 2 (n)\ <2C5 ( / ^-dr) . (2.22) 



r 2 .,2 \ h 
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However, in view of (2.18) and Fatou's lemma, we have 

-R r R 



/ ru 2 r dr < liminf / r(\u m } r ) 2 dr, (2.23) 
Jo m ^°° Jo 



R i , . .. . r R i 



r m->oo J r 



-w 2 dr < liminf / -w^dr. (2.24) 



In particular, in view of (2.18) again, we see that ^u 2 G L(0, R). Therefore the right-hand 
side of (2.22) tends to zero as r\,r 2 — > 0, which implies that the limit 

r] = \imu 2 (r) (2.25) 

r— s-0 

exists. Since ^u 2 G L(0, R), we must have r] Q = 0. Hence the boundary condition u(0) = 
is achieved. 

Summarizing the above results, we conclude that the function u obtained as the limit 
of the minimizing sequence {u m } for the problem (2.5) satisfies u(0) = u(R) = 0, u(r) > 
for all r G [0, R], E{u) < oo, and 

I(u) < liminf I(u m ), P{u) = lim P(u m ) = P . (2.26) 

Thus, u is a solution to (2.5). Consequently, there is some j3 G K. such that (w,/3) verify 
the boundary value problem (2.1)-(2.2). 

If there is a point r G (0, R) such that w(r ) = 0, then u r (r ) = since r is a 
minimum point for the function u(r). Applying the uniqueness theorem for the initial 
value problem of ordinary differential equations, we have u(r) = for all r G (0, R), which 
contradicts the fact P(u) = Pq > 0. This proves u(r) > for all r G (0, R). 

(ii) Let (u, j3) be a solution pair just obtained. We next study the quantity j3 in (2.1). 
As a preparation, we establish 

lim inf {ru(r) \ u r (r) \ } = 0. (2.27) 

Suppose otherwise that (2.27) is not valid. Then there are some eo > and r G (0, R] 
such that 

ru(r)\u r (r)\ > e , r e (0,r ), (2.28) 

which leads to 

oo = / — dr < u\u r \dr < I / -m 2 dr I I / r-u 2 dr ) , (2.29) 
Jo ?" io Vio r J \J J 

which contradicts with E(u) < oo. So (2.27) is valid. 

From (2.27), we can find a sequence {r m } so that r m — > as m — > oo and 

lim {r m u(r m )u r (r m )} = 0. (2.30) 
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Multiplying (2.1) by u, integrating over [r m ,R], letting m — >■ oo, and applying (2.30), 
we obtain r r r 2 

ft jf ra 2 dr = (rVu 2 + m 4 ) dr - ^!L w 2 + r ^ dr ^.31) 

Let m be any absolutely continuous function satisfying E(u ) < 00, the boundary 
condition (2.2), and -P(wo) = -fo- Since w solves (2.5), we have I(u) < I(u ). As a 
consequence of this observation, we have the bound 

J o (jT u2 + ru ^j dr - J q [rVu 2 + ^u 4 ^dr + 2I(u ). (2.32) 
Inserting (2.32) into (2.31), we obtain 

^-/3P > -21 (u ) + I [ R ru 4 dr. (2.33) 

To estimate the right-hand side of (2.33), we set R = 2a for convenience and define 

( -J, < r < a, 

Mr) = { a (2-34) 
[ J(2a-r), a < r < 2a. 

Therefore, after some calculation we have 

f 2a An 
P = 2tt ru 2 dr = —a 2 b 2 , (2.35) 
Jo 3 



2a 



r{u' Q {r)) 2 dr = 2b 2 , (2.36) 



f 1 

Jo r 

I 



dr = 26 2 (21n2-l), (2.37) 

r 

2a 



ru 4 dr = -a 2 b 4 . (2.38) 



Using (2.35)-(2.38) and (2.8), we get 

/(tio) < ^ R {rK(r)) 2 + ^-^}dr + i^rV- M 2 dr 

< b 2 (l + n 2 (2 ln2- I) + ^V - a 2 - -^a 2 b 2 J . (2.39) 



Besides, applying the Schwartz inequality, we have 

\ 2 R 2 f R 

rw 2 dr) < — / ru 4 dr. (2.40) 
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Thus, in view of (2.33), (2.39), and (2.40), we arrive at 



2tt 
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^f3P > 2b z I — a 2 b 2 - 



45 



— V ~R 2 + l+ n 2 (2\n2-l) 



(2.41) 



Inserting R = 2a and a b = 3Po/47r, we obtain the lower estimate for (3: 



12/7 



R 2 V60tt 



Pn 



— l/-p 2 + l + n 2 (21n2- 1) 



(2.42) 



(iii) We next derive the sufficient condition stated to ensure f3 < 0. 
In fact, inserting (2.13) into (2.31), we have 



2^° 5 



< 1 -"f™>-jf(T 



p2 

rV ) w 2 dr. 



r A-K 2 er 

For convenience, we may set £ = 1 in (2.43). Thus, whenever the inequality 

n 2 >^P 2 + r 2 V + (r), re[0,R], 



(2.43) 



(2.44) 



is fulfilled, we can conclude with /3 < since u(r) > for r £ (0, i?). 

(iv) We now consider nonexistence. For any admissible function u, we may view 
uas a radially symmetric function defined over M 2 with support contained in the disk 
Dr = {(x,y) £ M 2 | x 2 + y 2 < R 2 }. Hence, applying the classical Gagliardo-Nirenberg 
inequality over M 2 , we deduce 



R 



R 



ru dr < Att / ru dr / ru r dr, 



R 



(2.45) 



From (2.31) and inserting (2.45) with P(u) = P , we have 

= J (rVu 2 + ru 4 ) dr - J (^yU 2 + ru 2 . + firu 2 ^ dr 
< (2P -l)jT ru 2 r dr-J^ (j^ + ft - V^J ru 2 dr. 



(2.46) 



Therefore, when the conditions 

Pd S tt, -5 



V(r) + p>0, r£(0,P], 



(2.47) 



are imposed, we arrive at u = 0, as anticipated. Thus, in this situation, the problem 
consisting of (2.1)-(2.2) has no nontrivial solution. 
The proof of Theorem 2.1 is complete. 
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3 Vortices as saddle points 

In this section, we study the existence of optical vortices which are the solutions of the 
boundary value problem (2.1)-(2.2) as the saddle points of the action functional 

h{u) = l£ {ru 2 r + yu 2 + (/3- V(r))ru 2 - ^ dr, (3.1) 

with |n| > 1. We shall use a min-max theory approach. Suggested by the discussion of 
the previous section, we introduce the function space H which is the completion of the 
space X = {u G C^O, R] \ u(0) = u(R) = 0} (the set of differentiate functions over [0, R] 
which vanish at the two endpoints of the interval) equipped with the inner product 

(u,v)h — J (^ru r v r H — uv^j dr, u,v G H. (3.2) 

As seen in the discussion of the previous section, as a Hilbert space, H may be viewed as 
an embedded subspace of Wq ,2 (Dp>) consisting of radially symmetric functions such that 
any element u G H enjoys the desired property u(0) = 0. 

In order to simplify the presentation of the study here, we assume that (3 satisfies 

P > V + = max{\/ + (r) | r G [0, R]}. (3.3) 

Recall that a (^-functional / : H — > R is said to enjoy the Palais-Smale condition if for 
any sequence {u m } satisfying the properties (i) I(u m ) — > a as m — > oo, and (ii) I'(u m ) — > 
as m — > oo as a sequence in the dual space of H, one can extract a subsequence from 
{u m } which converges (strongly) in H . 

As an initial step, we have 

Lemma 3.1. The action functional I p given in (3.1) satisfies the Palais-Smale condition. 

Proof. It is straightforward to see that the functional (3.1) is C 1 over H . 
Let {u m } be a sequence in H satisfying the properties 

— >■ a, m — > oo, (3.4) 
{r\u m \ r v T + y« m f + C# - V{r-))ru m v - rw^r j dr 

< SmIMlif, £m > 0, V G H, (3.5) 

where e rn — > as m — > oo. In (3.5), we may take v = u m to get 

jT ru 4 m dr<J^ ^ri[u m ] r ) 2 + ^u 2 m + (P-V(r))rul^ dr + e m \\u m \\ H . (3.6) 
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\^m)(v)\ = 



On the other hand, in view of (3.4), we may assume without loss of generality that 
Ip{u m ) < ol + 1 for all m — 1, 2, • • •. Hence, applying (3.6), the assumption (3.3), and a 
simple interpolation inequality, we find 

1 f R ( n 2 \ 1 

2(a + l) > g i ( r (l ttm W 2 + 7 M ™ + ^ ~ ^W)™! I dr - 2 £m||Mm||ff 

> \\Wm\\ 2 H -\s 2 m , m = l,2,.... (3.7) 

In other words, {w m } is a bounded sequence in H. Without loss of generality, we may 
assume that {u m } weakly converges to an element u G H as m — > oo. It is clear that 
u m — > u as m — > oo strongly in any U'(Dr) or L p ((0, i?), rdr) (p > 1). Thus, letting 
m — ?■ oo in (3.5), we arrive at 



f R f n 2 ) 

/ < ru r v r H uv + (j3 — V(r))ruv — ru 3 v > dr = 0, Vi> G i?. 

io I r J 



(3.8) 



Let v = u m — -u in (3.5) and (3.8) and insert the resulting (3.8) into the resulting (3.5). 
We have 

J ^r([u m - u] r ) 2 + y(u m - uf + (p - V)r(u m - u) 2 - r(uf n - u 3 )(u m - «)| dr 

< £m\\u m - u\\ H . (3.9) 
As a consequence of (3.9) and [5 > Vq , we obtain 

\\u m -u\\ 2 H < e m \\u m -u\\ H + / r|(u^ - w 3 )(w m - u)\ dr, m = l,2,---, (3.10) 

which immediately implies that u rn — > u strongly in H as m — > oo, as desired. 

We next identify a mountain-pass structure through the following two lemmas. 
Lemma 3.2. There are constants K > and C > snc/i £/ia£ 

inf{/g(-u) | \\u\\ 2 H = K} > C . (3.11) 
Proof. For any constant K > 0, let u G H satisfy \\u\\ 2 H = K > 0. From (2.12), we have 

r R r R r R r R u 2 

/ rw 4 dr<4/ rdr / ru 2 Ar —dr<2R 2 K 2 . (3.12) 
</o -/o Jo Jo r 

Applying (3.12) in (3.1), we find 

Ip(u)>±(K-R 2 K 2 )=f(K). (3.13) 

However, the maximum of the function / in (3.13) is attained at K = ^ which gives 
us the value f(K ) = So, in conclusion, we have the lower bound 

«»)>^, II* = 5^. < 3 - 14 > 

which establishes (3.11). 
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Lemma 3.3. For any constant K > 0, there is an element v G H satisfying \\v\\ 2 H > K 
and Ip(v) < 0. 

Proof. With R = 2a, we will see that we can use the function u defined in (2.34) as a 
test function. 

For this purpose, we first show that u G H . To see this, we need to prove that Uq can 
be obtained in the limit from a sequence of functions in X under the norm of H. In fact, 
for any < e < a, we can define 

(uo(r), r G [0, a — e) U (a + e, 2a], 
(3.15) 
Q £ (r), r e[a-e,a + e], 

where Q £ (r) is taken to be a quadratic function satisfying 

Q £ (a ± e) = u (a ± e), Q' e (a ± e) — u' Q (a ± e) . (3.16) 

Matching these conditions, we find 

Q £ (r) = - 7 ^(r 2 -2ar+ [a-e] 2 ), re[a-e,a + e], (3.17) 
2ae 



which enjoys the bounds 

b, , - „ , , .h 
a 

Consequently, from (3.18) we have 



(a-e)<Q £ (r)<-(a- 6 -), \Q' £ {r)\ < -, re[a-e,a + e]. (3.18) 



lim^ a+£ (r[Q' £ {r)] 2 + ig e 2 (r)^) dr = 0. (3.19) 

Therefore {u £ } is a Cauchy sequence in H as e — > whose limit is clearly uq in view of 
the definition of {u £ } given in (3.15). This proves u G H. 
Using the results (2.35)-(2.38), we have 

Htxoll^ = 46 2 ln2, (3.20) 
Wo) < ^jf fl {r([«o]r) 2 + y«g + (/3 + V-)r«g-^}dr 

< 6 2 ^l + n 2 (21n2-l) + ia 2 (/3 + \/ -)-^a 2 6^. (3.21) 

From (3.20) and (3.21), we see that for any K > we may choose b > sufficiently large 
to get Ip(u ) < and ||uo||h > K- 
Thus the lemma follows. 

It is interesting to note that (3.21) implies that Ip(u ) — > — oo as b — > oo. In other 
words, the action functional (3.1) is not bounded from below over H which prevents a 
direct minimization approach to the problem. Indeed, we are now prepared to obtain 
a nontrivial solution of the boundary value problem (2.1)-(2.2) as a saddle point of the 
functional (3.1) in the following theorem. 
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Theorem 3.4. For any (3 > V + , \n\ > 1, and R > 0, the problem (2.1)-(2.2) has a 
nontrivial solution over the interval [0,R]. Moreover, such a solution may be obtained 
from a min-max approach applied to the action functional (3.1). 

Proof. Let Ip be the action functional (3.1). Then Lemma 3.1 says that Ip satisfies the 
Palais-Smale condition. Let K,Cq > be the constant stated in Lemma 3.2. Using 
Lemma 3.3, we can find some uq G H such that ||uo||/f > K and Ip(uo) < 0. Denote by 
T the set of all continuous paths in H that link the zero element of H to u : 

T = {ge C([0, 1]; H) \ g(0) = 0, <?(1) = u } . (3.22) 

Therefore there is some point t g G (0,1) such that ||<7(i ff )||# = K. From the classical 
mountain-pass theorem (cf. Evans [12]), we know that 

c = inf max Ip(g(t)) > C , (3.23) 

is a critical value of I p. In other words, there is an element u G H satisfying Ip{u) = c 
which is a critical point of Ip. Of course, u is nontrivial. That is, u cannot be the zero 
element of H. 

Recall that in Theorem 2.1 the condition (2.10) is obtained to ensure the wave prop- 
agation constant f3 to assume a negative value. Although (2.9) states a lower estimate 
for (3, no condition has been obtained to ensure /3 > for the solution of the constrained 
minimization problem (2.5). Theorem 3.4, however, complements Theorem 2.1 in that it 
gives us a family of nontrivial solutions realizing arbitrarily prescribed parameter in the 
entire interval |V + ,oo) for any vortex charge |n| > 1 and R > 0. That is, our existence 
result indicates that j3 may take any positive value above or at V + . 



4 The defocusing case when 5 = — 1 

If we have s = — 1 in (1.1) instead, then the action functional (2.3) is replaced by 

1 f R ( n 2 r *) 

I(u) = -J lru 2 r . + — u 2 -rV(r)u 2 + -u 4 j dr, (4.1) 

the difficulty with the quartic term, which was negative before, disappears, and the con- 
strained minimization problem (2.5) is easily solved, which gives us a solution to the 
associated equation 



n 2 

(ru r ) r u + r(V — u 2 )u = (3ru, (4.2) 



for some /3 G R. As before, this equation leads us to the relation 

f R f R ( n 2 _ \ 

J ((3 — V + )ru 2 dr = — J \r-u 2 . H u 2 + rV u 2 + ru 4 J dr. 



(4.3) 
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Thus, if j3 satisfies (3 > V + , then u = 0. Therefore, regardless of the value of R, the 
problem prevents the existence of a nontrivial solution for sufficiently large propagation 
constant (3. This conclusion is in sharp contrast to that in the case when s = +1 stated 
in Theorem 3.4. 

In general, the simple relation (4.3) clearly indicates that it is natural for (3 to take 
negative rather than positive values. For example, using the Poincare inequality over D R , 



R 2J (R 

ru dr < — 



R 



ru 2 dr, 



(4.4) 



where r (~ 2.404825) is the first positive zero of the Bessel function J , and (2.20), we 
obtain from (4.3) the result 



R 



ru 2 dr < V + 



< - 



R nR 

ru 2 dr — 
o Jo 



n 



ru r 



u + rV u l 



ru 



dr 



Consequently, we obtain 



r_o_ 
R 



- V + ^j J ru 2 dr - J ; 



ru dr. 



n 



R 2 



(4.5) 



(4.6) 



In particular, we have (3 — > —oo as \n\ — > oo. 

An example of the defocusing case s = — 1 is the study carried out in [18] where V takes 
the form V = pJf(br) given in terms of the Bessel function J\ and positive parameters 
p,b. Thus V~ = 0. Another example of the case s = — 1 is in the lines of the studies 
[24, 30]. There, although V is not radially symmetric, it is non-positive valued, V + = 0. 
Here, assuming V is radial as well as non-positive, then (4.3) indicates that (3 < is the 
only possibility. 

Note that the (z-direction) angular momentum of the obtained stationary vortex wave 
in view of (1.2) has the simple but elegant expression 



Im / (ijfdeijf) rdrd# = 2-kh J ru 2 dr 



nP. 



(4.7) 



u(r)e , where u solves (1.3) and satisfies u(0) = 0, and A 



Write 0(x, y) 
Then satisfies 

A0+(\/ + S |0| 2 )0 = /30, (4.8) 

away from the origin of M 2 . The condition u(0) = ensures that the origin is a removable 
singularity [35] such that when V is an analytic function of (x,y), so is 0. Consequently, 
in this situation u vanishes at r = like r n for an n-vortex solution as in the classical 
Ginzburg-Landau equation case [8, 14, 25]. 

For the focusing case s = +1 with a non-positive potential (cf. [34]), we have V + = 
and the statements of our results simplify considerably. For example, for the solution pair 
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(u,j3) obtained in Theorem 2.1 to have the property j3 < 0, it suffices that the vortex 
number n satisfies the condition 

N > £. (4.9) 
Moreover, applying Theorem 2.1 (iv), we see that there is no nontrivial solution satisfying 

P(u) < \ and n 2 > -r 2 /3, re[0,R}. (4.10) 

In particular, we conclude that there is no nontrivial solution with P{u) < \ and j3 > 
0. Besides, in this case Theorem 3.4 becomes an existence theory for any prescribed 
propagation constant j3 > 0. 
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